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Abstract We calculate zero temperature Green's function, the density- 
density correlation and expectation values of a one-dimensional quantum 
particle which interacts with a Fermi-sea via a (5-potcntial. The eigenfunc- 
tions of the Bethe-Ansatz solvable model can be expressed as a determinant. 
This allows us to obtain a compact expression for the Green's function of the 
extra particle. In the hardcore limit the resulting expression can be analysed 
further using Painleve V transcendents. It is found that depending on the 
extra particles momentum its Green's function undergoes a transition of that 
for hardcore Bosons to that of free Fermions. 
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£*"»« . 1 Introduction 

Many intriguing questions are related to the dynamics of one-dimensional 
systems QQQ]- Landau's Fermi-liquid theory fails for Fermionic systems in 
one-dimension and the Luttinger- liquid model applies [5] . The free Fermi-gas 
as a fixed point of the renormalization group is unstable in favor of collective 
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Bosonic excitations. This makes it difficult for an injected Fermionic particle 
to deposit its excess energy and yields largely different electronic transport 
properties in one-dimension [JJ. Corrections to the Luttinger-liquid model 
were studied recently experimentally [3J and theoretically [2J. On the other 
hand the unusual slow decay of a one-dimensional interacting Boson system 
has been observed [5] and was related to the exact classical and quantum 
solvability of the Licb-Liniger model [HIE]. This model was studied using 
Luttinger liquid theory in |S|. 

If in a Id gas different species of particles are present the theoretical de- 
scription becomes more difficult, since the many-body wave function trans- 
forms according to a higher dimensional representation of the permutation 
group. This results in additional quantum numbers apart from the asymp- 
totically free momenta, which have no classical analogue. Although for 5- 
interaction the eigenfunctions have been constructed via the nested Bethe- 
Ansatz [MTU] the problem of calculating Green's functions, respectively cor- 
relation functions are only solved in some special cases [TTllT2l[T3] . In spin-| 
systems the number of additional quantum numbers is infinite in the thermo- 
dynamic limit. However, if one spin polarization, say spin up, is carried by 
just one particle, only one additional quantum number appears. As we will 
show it can be identified with the free momentum of the extra particle. The 
eigenfunctions of this model as well as the extra particle's density-density 
correlator were calculated by McGuire jl41IT5] . He focussed on the ground 
state, i. e. the additional quantum number is not present in McGuirc's work. 
If one is interested in the dynamical properties of the distinguishable particle 
it is crucial to include it in order to describe excitations of the extra particle. 

We consider McGuire's model where the distinguishable particle is al- 
lowed to carry an arbitrary momentum and study its response to the in- 
teraction with a Fermi-sea, which is assumed to be at zero temperature. It 
might be viewed upon as an extremely imbalanced mixtures of two particle 
species. Their study was pioneered by Lai and Yang [TB] and has revived 
boisterously in the recent years both theoretically [1511171 118,19,20] and ex- 
perimentally [2T1I22U23] . 

We write the full Bethe-Ansatz many-body wave function as a deter- 
minant and are thus able to derive compact determinantal expressions for 
the extra particle's equal time Green's function (or likewise reduced den- 
sity matrix). In the Tonks-Girardeau regime of infinite interaction strength 
we employ the powerful methods of Tbplitz determinants and Painlcvc equa- 
tions, which were applied earlier to hardcore Bosons [23 > 25 , 26, 27, 28] . We find 
that in the hardcore limit the particle's expectation values and its Green's 
function depend crucially on its momentum. If the momentum is right at the 
Fermi edge the extra particle behaves just like an additional Fermion of the 
sea. The energy shift is zero and its single particle Green's function is that of 
a free Fermion. If the particle's momentum is in the core of the Fermi-sea, 
the energy shift is finite and the Green's function is identical with that of a 
hardcore Boson. As the particle's momentum varies from the Fermi edge to 
the core, the particle undergoes a transition from a free Fermion to a hard- 
core Boson, manifested in an algebraic asymptotic decay of the real space 
Green's function lim^i G(x) ~ with an exponent j3 changing from one 
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to 1/2. A small distance expansion of the Green's function and numerical 
work show that for finite interaction strength this picture does not change 
qualitatively. 

The article is organized as follows: In Sec. [2] we introduce the model. 
Due to its importance in the present context, we give a detailed discussion 
of the additional quantum number appearing in the Bcthc-Ansatz equations. 
Furthermore expectation values are calculated. In Sec. [3] we derive a determi- 
nantal expression for the Green's function and analyze it further in the limit 
of small distances and in the hardcore limit. In Sec. 2] the density-density 
correlation function is derived. We summarize and conclude in Sec. [5] 



2 Model 

We consider a N + 1 particle system in one dimension where N particles are 
identical spinless Fermions to which we refer as environment or Fermi-sea. In 
addition we assume the presence of a single particle which is distinguishable 
from the former ones. The particles interact via a repulsive 5-potential. Due 
to the Pauli-principle, it acts only between the extra particle and the particles 
of the Fermi-sea. The particle masses are all equal. Choosing the units such 
that h — 1 and m = 1/2 the Hamiltonian reads 

N o2 Q2 n 

Here the coordinates Xj, j = 1, . . . , N denote the positions of the N identical 
Fermions while y refers to the position of the distinguishable particle. Fur- 
thermore c > denotes the interaction strength and the factor 4 is included 
for convenience. The model is exactly solvable via Bethe's-Ansatz rTllll5| . 
However the original form of the wave function constructed by this method 
is inconvenient for explicit calculations. 



2.1 Eigenfunctions and Bcthc-Ansatz Equations 



Denoting an cigenfunction of the Hamiltonian ((T|) by W(xi, . . . , xjy, y), we 
require it to be antisymmetric in the first ./V arguments. This and the fact 
that there is no interaction between the particles of the Fermi-sea suggest to 
write ^(xi, . . . , xn, y) as a determinant 



W{x u ...,x N ,y) = C w det Aj(x t - y)e tk ^' 
Aj{x) = i(kj — A) + csgn(x). 



j = l....,N + l 

l=l,...,N 



(2) 



Note that \&(x\, . . . ,xn, y) has all required symmetries i.e. is antisymmet- 
ric under the exchange of all quasimomenta, antisymmetric in the particle 
postions xi , I =, 1, . . . , N and has no well defined symmetry when the distin- 
guishable particle and a particle of the Fermi-sea are exchanged. As shown in 
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Sec. 13.21 the normalization constant Cm is for periodic boundary conditions 
determined by (see Eq. ([73])) 



N+1 N+1 



\C N \- 2 = NIL n [Lm - Af + c> + 2C/L)] g L{{kj _ A)2 + c2 + 2c/ L y 



(3) 



In Eq. ^ {kj}j = \.... t N+i are the " quasimomenta" and A is an additional 
parameter (not present in [HUTS] ). 

The wave function © is an eigenfunction of the Hamiltonian (JlJ to the 
eigenvalue 

JV+l 

E =Y. k r ( 4 ) 

3=1 

Furthermore the eigenvalue to the center of mass momentum operator is 

N+1 

K=Y d k j . (5) 

3=1 

The form ((2J) is crucial for our approach, since it allows us to employ powerful 
methods of matrix algebra to manipulate determinants. Equations ©-([S]) can 
be proved straightforwardly by acting with H on W. 

In order to acquire results for finite particle density we impose periodic 
boundary conditions 

<?(£i, . . . ,xi, . . . ,x N ,y) = &(x!,...,xt + L, . . . ,x N , y) , I = 1, . . . , N, (6) 
>P(x 1 ,...,x N ,y) = &(x lt ...,x N ,y + L) . (7) 

The first condition yields the transcendental equations 



e 



tki L 



i(kj — A) + csgn(ai — y) 



, j = 1,..., JV + l (8) 



i(kj — A) + csgn(x + L — y) 
which after taking the logarithm on both sides translates to 

( k — A\ 

kjL = 2nrij - 2 arctan f J , j = 1, . . . , N + 1 . (9) 

As can be seen by considering the limit c — > of Eq. (|9|) the numbers n 3 - 
have to be half- integers and — n/2 < arctan(x) < n/2. The second condition 
yields 

N+1 _ 

kjL = lithj + 2 ^2 arctan ( — J, j = 1,...,N+ 1, (10) 

where fij are half-integers for N odd and integers for N even. In the following 
we will for convenience always assume N to be odd. 
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Combining Eq. ([8]) and Eq. (|10[) yields quantization rules for the quasi- 
momenta kj and A which are known as Bethe-Ansatz equations 

kjL = 2tt (^rrij + - 2 arctan ^ , j = 1, . . . , N + 1, (11) 

N+l f k - A\ 
2ttJ a = 2 ^ arctan ( J . (12) 

Now the nij are unequal integers. According to Eq. (|12p the integer Jyi is 
bounded by 

- ^ s * < ^. 

The quantities , are the quantum numbers labeling an eigenstate of the 
system. For J a = A = an expression equivalent to Eq. (jllj) was derived in 
Ref. [13]. For a derivation of Eqs. ([TO and (JT2J) from the Bethe-Ansatz wave 
function we refer to Ref. [7]. 



2.2 Solution of Bethe-Ansatz Equations and Ground State 



The Bethe-Ansatz equations are coupled algebraic equations which in general 
do not permit a closed solution. In our case the situation is more favorable, 
since the equations for kj only couple via A. For c — ¥ and for c — > oo they 
can be solved exactly. 

Let {jtIj}^! jv+i be a set of ordered integers nij > m.j + i. Then in the 



limit c 

by 



+ a solution to the Bethe-Ansatz equations (TTTj) and (|T2"j) is given 



2tt , 
2tt 



jV+1 

2 i 



(14) 



— K- + i) ,j>j^ + iy±i 



and A € [m^m^i + 1]. The superscript ± indicates that kj is approaching 
its value from above/below. For c-)0 + the solution for A is in general not 
unique. We come back to this point below. 

Since the quantum numbers nij and J a for the ground state are the same 
for c — and c > 0, the quantum numbers nij of the interacting ground 
state can be extracted from Eq. (Till) 



{ m j}j= 



1,...,JV+1 



fJV-1 



N+ 1 



(15) 



This yields for c = the set of quasimomenta 

27 



{%}. 



<jjj=l,...,JV+l 



2?r r 



o + ,o-,±i, 



,±- 



N-l 



(16) 
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Since the distinguishable particle can occupy the same state as a sea-particle, 
in the ground state the single particle quasimomentum with the lowest energy 
is double occupied. Thus for the ground state J a = 0. The allowed interval 
for A shrinks to a point and the solution for A becomes unique A = 0. The 
same happens for any other state for which rnj A+(N+1)/2 and m (j A+(N+1)/2 +i) 
are adjacent integers. 

In the hardcore limit c — > oo Eqs. (fTTj) and (JT2J) become 



where the dimensionless factor A takes into account that A scales like c for 
c — > oo. From the definition in Eq. (|18[) it is seen that A can range between 
— oo and oo when J a varies according to Eq. (|13[) . 

For finite c the Bethe-Ansatz equations can be calculated numerically. As 
c increases the quasimomenta kj evolve smoothly from 2iv(rnj +6(A — rrij))/L 
to 27!"(TOj + iSgn(jl— rrij))L, i. e. from multiple integer to multiple half-integer 
values of 2ir/L. The additional parameter A evolves smoothly from its initial 
value to ±oo as the interaction increases. 

In the following we assume that the sea-particles are in the ground state 
but the extra particle may occupy an arbitrary single particle state. Two 
cases have to be distinguished: 

Either the quasimomentum of the extra particle lies inside the Fermi-sea. 
Then the quantum numbers nij of the state are given by the set of integers 
in Eq. (|15p. The non-interacting extra particle shares its momentum with 
one of the particles of the Fermi-sea. The additional quantum number J a 
indicates the double occupied quasimomentum. For J a = (N — l)/2 or J a 
= —(N— l)/2 the Fermi-momentum of the sea is double occupied. The non- 
interacting eigenfunctions for a given set of quasimomenta {fcj}j=i,....Ar+i is 
uniquely determined and the energy-eigenvalue is non-degenerate. Thus for 
c = the energy increases quadratically with J a- 

In the second case the quasimomentum of the extra particle lies outside 
the Fermi-sea. The quantum number J a takes its value at the upper edge 
Ja = (N+l)/2 or at the lower edge Ja = — (N+l)/2 of its spectrum. For van- 
ishing interaction all N + 1 quasimomenta are different. If the highest /lowest 
quantum number is mi = (N — l)/2 respectively mjy+i = —(N + l)/2 the 
extra particle's momentum borders the Fermi-sea. If mi > (N — l)/2 or 
tojv + i < — (N + l)/2 the extra particle's momentum is outside the Fermi 
sea. For a given set of quasimomenta there exist N + 1 orthogonal eigenfunc- 
tions since the extra particle might carry any of the fc n 's without changing 
the systems energy. These eigenfunctions are distinguished by J a- Thus at 
the Fermi-edge the energy becomes independent of Ja and the density of 
states has a singularity. The different cases are illustrated in Fig. [T] 




(17) 



(18) 
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I) J a = 0, mi = 2 



II) J A = -3, mi = 2 



- » i ♦ t ■ 

3 2 1 0... 
Ill) J A = -2, mi = 2 
t 

0~vJ • • • • • H — 

3 2 1 0... 







3 2 10. 



IV) J A = -3, mi = 4 



3 2 10. 



Fig. 1 Sketch of the occupied non-interacting one-particle states for N = 5. A 
black circle denotes a state occupied by a sea-Fermion, an empty circle denotes 
an empty state. The grey circle denotes the state, which is occupied by the extra 
particle. Case I) is the ground s tate configuration. In case II) and III) the quantum 
numbers rrij are given by Eq. I|15p . In case II) the extra particle's momentum is 
just outside the Fermi-sea in case III) it is just inside. 



Summation over Eq. (|lip from j = 1 to j = N + 1 yields in combination 
with Eq. (fT2j) for the overall momentum of the system 



JV+1 



/N+l 



Ja 



(19) 



If the set {m,j}j = 



j}j=l,...,N+l 



is given by Eq. (|15|) the sum vanishes and the 
quantum number Ja — is directly related with the center of mass 

momentum K. The latter in turn is identified with the momentum of the 
extra particle in the lab frame where the Fermi-sea is at rest. 



2.3 Thermodynamic Limit 

For a dense set of quantum numbers as in Eq. (|15p the density of states 
g(k) = dm(k)/dk can be derived by taking the derivative of Eq. ([TT]) 

Q( k ) = ^~ + -T, 7- ( 20 ) 

w 2tt 7t (k- A) 2 + c 2 y ' 

This expression has a rather natural interpretation: The first term corre- 
sponds to the density of states of a non-interacting Fermi-sea. The second 
term appears due to the presence of the distinguishable particle. It has the 
form of a Lorentzian distribution centered around A with width c. For c — > 
it becomes a (5-function at k = A. 

In the thermodynamic limit, where N, L — > oo such that the particle 
density p = N/L remains finite there are different choices for the scaling of 
the interaction strength 

cL- 1 

c(ttN/L) . (21) 

oo (harcorc limit) 
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The Lorentzian part in Eq. ((20]) is outscaled by the part due to the Fermi-sea 
in all three cases. This implies that in the thermodynamic limit the single 
particle does not affect the spectrum of the Fermi-sea. The effect of the 
particle onto the bath should therefore be negligible. This is of course an 
essential requirement for a bona fide dissipativc quantum system. 
The system's total energy is 



+k ¥ 



E 



dkk 2 g{k) . 



(22) 



For the energy shift due to the interaction this yields up to corrections of 
order 1/N 



E — En 



kp 



kp 



A 2 



-v(kp, c, K) 



2tt 
~L 



J A 



(23) 



2k F + Aln 



where the Fermi-momentum is kp 

7T 2 N 3 



+ (A- kp) 2 
+ (A + k F ) 2 

tt P + O(N ) and 



En — 



3L 2 



27T 

T Ja 



O{N ) 



(24) 



corresponds to the energy of the non-interacting system. In the thermody- 
namic limit where the quasimomcnta are distributed according to Eq. (TT51 . 
it proves useful to define the quantity K = K/kp = 2Ja/N. It corresponds 
to the overall momentum in units of kp and measures the position, where 
the extra particle's momentum is located in the Fermi-sea . For K = the 
extra particle's momentum is in the center, for K = ±1 it is at the upper or 
lower edge. 

Since it appears frequently troughout the following we introduced in 
Eq. (|2"3"|) the abbreviation 



v{kp,c, K) = 



arctan 



A 



arctan 



A 



(25) 



Note that the parameter A in Eq. (|2"3"]l is not independent but relates via the 
second Bethe-Ansatz equation ([12"]) to J a, kp and c 



J A 



— / dkg(k) arctan 

7T J 



N 

— / ax arctan 
2tt 



k-A 



A/kf 



From this equation follows in particular that A\ c= 



O(N ) . 
= -2nJ A /L. 



(26) 
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Although the integral can be evaluated the resulting transcendental equa- 
tion can not be solved analytically for A. Therefore the energy shift can in 
general not be expressed as function of the three independent quantities J a, 
k F and c and has to be calculated numerically. For the ground state where 
J A = A = 0, the situation is more favorable and the expression (l23l) simplifies 
to 



E-En = 



ckp + kp arctan 



— c arctan 



k\ 



The mean interaction energy 

(V) = (v \4cJ2Hxn-y)\&o) 



N 



n=l 



can be obtained using Pauli's trick (see Ref. [53]) as 



(V)=c 



dE 4c 

dc 7T 



kp — c arctan 



(27) 



(28) 



(29) 



The interaction energy is given by the interaction strength times the prob- 
ability to find a sea-particle at the same position as the extra particle. The 
latter is identical with the equal time density-density correlation function 
R(x, y) at x — y — 0. Thus we have, by definition, 



(V) = 4cLR(0,0) . 
Equating Eqs. J2SJ) and ([201) yields 



#(0,0) = 



iV 

p 



1 — c arctan 



(30) 



(31) 



Thus the probability to find a sea particle at the position of the extra particle 
vanishes as c -2 for c — > oo. Therefore the interaction energy vanishes in the 
weak coupling limit as well as in the strong coupling limit. 

In Sec. [31 where the full density-density correlation function is studied, 
we derive an expression for R(0, 0) for J a ^ 0, see Eq. ([90]) . This allows us to 
calculate the interaction energy in this case as well. In Fig. [5]the ground state 
energy shift and interaction energy are plotted versus the interaction strength 
for three different values of K. It is seen that the interaction energy has a 
maximum as predicted. The ground state energy shift increases monotonously 
and saturates at a value E meLX = k F — K 2 . This means that for the extra 
particles momentum located right at the edge of the Fermi-sea the energy 
shift vanishes for arbitrary interaction strength. 

Using the wave function the mean expectation value of the extra 
particles energy can be calculated. Up to corrections of order 1/N we obtain 



(Pi) 



A 2 



v~ 1 (k F ,c,K) (2fc F c + cylln 



c 2 (k F 



Af 



c 2 (k F + Ay 



(32) 
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Fig. 2 Interaction energy (full lines) and energy shift (dot-dash lines) as function 

of c. The values of K are K = (black), K = 0.5 (online red) and K = 0.9 (online 
blue). 



For the ground state this simplifies to 

«*> = *Ki^(Fi)- a )- (33) 

From Eq. (|33p it is seen that the extra particles energy increases monotoni- 
cally from (p^) = to (p^) = fc|/3 as c varies from zero to +oo. 



3 Equal Time Green's Function 



The equal time single particle Green's function with respect to a wave func- 
tion <P is defined as the iV-fold integral 



G(y,y') 



dxi 



(34) 



Likewise it can be interpreted as the reduced density matrix of the ex- 
tra particle. We stick to the notation equal time Green's function or sim- 
ply Green's function. Using the representation ([2]) for the wave function 
closed expressions can be obtained for G{y,y'). Expanding the determinant 
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in Eq. ([2]) with respect to the last column and using the normalization con- 
dition G(0,0) = 1/L we obtain 



N+l N+l 



G{y,y') = \C N \ 2 Y. EC" 1 



n+m e i(k n y-k m y') 



(35) 



n— 1 m=l 
L L 



/ dxi . . . / dx N dct [A^ajj - y)e %k * Xi } j= i n+i^t, 

J J 1=1,..., N 



det [A*(xi - y')e lk > Xl ] j=i,..., N +i& 

1=1,.. .,N 



Employing the general result 

J dxi . .. J dx N det [/j(^)].,-,;=i,...,jv det [hj{xi)] jl=1 _^ N (36) 



= AM dct 



dxhj(x)fi(x) 



j.l=l,....N 



one determinant can be replaced by its diagonal part and the ^-integrations 
can be performed. The resulting expression acquires the form 



N+l N+l 



G(y,y') = \C N \ 2 N\Y, £(-l) n+m e 



t(k n +k m )(,y-y')/2 



x det [Kji(y,y')] j=lt .... N+l9in , (37) 

l=l,...,N+l^tm 

where the matrix entries are given by 

L 

Kjl(y,y') = e-< k *- klKy+y ' )/2 [ dxAj(x - y)A* l (x - y')e l ^- kl)x . (38) 



The evaluation of Kji(y ) is straightforward but tedious. It yields 

K jt (y-) = [L((k 3 - A) 2 + c 2 ) + 2ct(kj - A)y~ - 2c 2 y~)} 5 fl (39) 
'(kj - k t )y- 



2c < cos 



h + kj -2A + z2csgn(y ) . / (kj - h)y 



kj ki 



sm 



(1-%) 



with y~ = (y — y'). The expression (j39|) makes explicit that G(y,y') is a 
function of the difference y~ only, as expected by translation invariance. 
After some further matrix algebra the r. h. s. of Eq. (|37[) can be can be 
expressed in terms of a single N x N determinant 



G(y,y') 



-Gi(y-), 



(40) 
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where we have introduced the interaction part G\ (y ) of the Green's function 
Gi(y~) = \C N \ 2 iVILdet [g nm - g nm +i - 9n+im + SW+im+i] „ im= i.... jA r( 41 ) 
and the quantities 

g nm = e-< k - +k ^~' 2 K nm {y-). (42) 

According to Eq. (|4"0")) we have G\(0) = 1. Separating the overall momentum 
and G\(y~) as in Eq. (|4*fJ)) is useful when exited states with K ^ are 
studied. The representation of the Green's function as a determinant is most 
convenient for a further numerical analysis. 

In Figs. [3] and 0] the real part of G\(y~) is plotted for a particle number 
of N — 15. The quantum numbers are chosen according to Eq. (|15[) and 
c = cL. Figure [3] shows the transition of the interaction part of the Green's 
function for the ground state as c varies from c = to c = oo, and Fig. 
2] shows Gi(y~) for exited states with J a ^ and fixed c = 50. Whereas 
in the ground state the curve is smooth and decays monotonously up to 
y~ /L = 0.5, for excited states wiggles develop and the function seems to 
become oscillatory. However, only for the highest value J a = 8, i. e. for the 
particle just outside the Fermi-sea, the curve has nodes. 



J A = A = 0,N= 15 
c = 

c = 20 /'" 

"c = 50 



c = ooy' 



0.2 0.4 0.6 0.8 1 

y~/L 

Fig. 3 Real part of Gi(y~) as function y~ /L for the ground state i.e. J a = and 
iV = 15. The values of c are c = (full line, black), c = 20 (dashed line, online 
red), c = 50 (dotted line, online blue) and c = oo (dot-dashed line, online green). 

In the following, we analyze the representation (|4"Tj) of Gj(y~) further in 
the hardcore limit. We show that the features, shown in Figs. [3] and |4] and 
described above persist in the thermodynamic limit. 



OA 



0.8 



0.6 



0.4 



0.2 
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c = 50,N = 15 




0.2 0.4 0.6 0.8 1 



y~/L 

Fig. 4 Real part of Gi(y~) as function y~ /L for c = 50 and JV = 15. The values 
of J a are J a = 4 (dot-dashed line, black), Xi = 6 (dotted line, online red), and 
J a = 8 (full line, online blue). 



3.1 Hardcore Limit 

For c y oo the wave function ([5]) becomes 



&{xi, . . .,XN,y) oc ( - + sgn(x i - y)j dot 



ikjXi 



e 



j = l,...,7V+l ' 
i=l,...,N 



(43) 



where we recall A = yl/c. For A = the expression resembles the wave 
function of N + 1 hardcore Bosons The only difference consists in the 
product of sign-functions: The wave function for hardcore Bosons is a Slater 
determinant multiplied not with a single product of sign-functions but with 
a double product n»<j s S n ( x i ~ x j)- For the equal time Green's function this 
difference is irrelevant and it is identical with those of hardcore Bosons. A 
large body of literature has been devoted to the studies of the latter [30ll24l 

EDEllllklllHlEllMllSa . 

On the other hand for |A| — > oo the wave function and therefore also 
the Green's function coincides with that of free Fermions. This is already 
reflected in the Bethe-Ansatz equations (|17[) and (|T8|) . The quasimomenta 
are integer multiples of 2n/L with an offset arctan(A)/7r which ranges from 
zero to ±i for |A| — > oo. 

Hence varying A from A = to A = oo we expect the extra particle to 
undergo a transition from a hardcore Boson to a free Fermion. This transition 
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should be reflected in the equal time Green's function. As there are two 
fundamental formulations of the Green's function for impenetrable Bosons 
we will go for two different ways to study it. The first approach is based on 
the theory of Toplitz determinants i30,24,251[2"T] and will be reviewed in the 
appendix [21 

The second approach is the description of the Green's function via solu- 
tions of a Painleve equation [25128 . In the following we are going to show 
that the above mentioned transition is for zero temperature described by a 
change in the initial condition of the solution of one and the same Painleve 
equation. 

We plug Eq. (f4"3"|) into the definition ([M)l of the equal time Green's func- 
tion. This yields 

N 

G(y, y') = \C N \ 2 J d[x] JJ(^ + sgn(ar, - y))(-i\ + sgn( Xl - y')) 



i=i 



x det 



j=l 7V+1 

l=l....,N 



det 



-ikjXi 



-ikj-y' 



3 = 1 N+ 

1=1. ...,N 



1 -(44) 



The crucial point is that for zero temperature, or - more precisely - if the 
Fermi-sea is in the ground state, the integral can be interpreted as an aver- 
age over the unitary group. We rewrite the last line of the equation above 
according to 



det 



j=l,...,N+ 
1=1,.. .,N 



1 det ^ 



— ikjXi \„ — ikjy' 



\ — arctan(A)?/ 



n - 2 

l<j<l<N 



j=l,...,N+l 
1=1,. ..,N 

irjxj - xi) 



(45) 



L 



N 

n 

i=i 



Sill 



( ttxi _ wy \ 
\ L L) 



sin 



( 



wy_ 
L 



In the present case the first product can be identified with the measure of 
the unitary group U(N) [5^ 



d/J.(U) oc sm 

l<j<l<N 



N 

n 

3=1 



xi)\ yj 2-irdxj 



(46) 



The Green's function can again be written like in Eq. (|40[) . where the part 
Gi, which includes the effects of interaction, is identified with the average of 
the function 

N 

ni 1 -^^) 2 ^^)^^'^) (l-e 2 ^^-'^)) (47) 
i=i 

with respect to the measure d[i(U). Here we have introduced the indicator 
function 



y (0 
A [y',y) 



1 for xi G [y',y) 
else 



(48) 
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and the parameter 



Averages over the unitary group of functions of the type (|4"7|) were studied by 
Forrester, Frankel, Garoni and Witte [371I271I28] and they have been related 
to solutions of the Painleve VI non-linear differential equation. Defining 

u = cxp — i — - — , (50) 



we follow Proposition 3 and Corollary 1 of Ref. [25] in order to deduce the 
relation 

u(u- l)-^\nGi(y~) _^_ /r =a N+ i(u) (51) 



i2ir B - /L 



for the interaction part of the Green's function. The function o\/v(it) is a 
solution of the particular Pyi er-form due to Okamato [38] 

- u 2 (u - 1) 2 K0 2 - (a N - (u - l)a' N + 1) (52) 
x [4a' N (a N -ua' N )- (N 2 - l)(a N -(u- l)a' N )] . 

The initial conditions arc fixed by expanding Gi(y~) for small distances like 
in Ref. [35] . We find that in terms of the variable u the solution of Eq. (|52l) , 
quested for, has the small distance expansion 

N 2 — 1 

a N {u) = ^-{u-l) 2 (53) 



(N 2 -l) (iN/(l + iX)-n) f 
•- (u - 1 



24. 

In the thermodynamic limit the expressions simplify further. First, Painleve 
VI converts in the limit N — > oo to the Jimbo-Miwa-Okamoto form [26,39 
of the Painleve V differential equation 

^W) 2 = (ov - xt/y + 1) {{a' v f -a v + xa' v ) , (54) 
where x = k-py~ ■ Second, the initial condition for <Jy{x) becomes 

oy(x) ~ -Y + lhx^ + — X ^°- (55) 
Thus the Green's function reads in the thermodynamic limit 

p iKx 

G(x) = —j—Gi{x) 

G I {x)=cxp(ixK+ [ dx'^p-\ . (56) 
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The long distance behavior of Gj(x) can be related to the short distance 
behavior (|55p . i. c. to A with a connection formula for Painleve equations. 
The case of Painleve V has been solved by McCoy and Tang [10114 11142] . We 
employ the results obtained in Ref. [H] for hardcore Bosons: 

Let the solution of Painleve V in the Jimbo-Miwa-Okamoto sigma form 
as given in Eq. (|54j) be regular at x = with the expansion 



a v (x) ~ -- + - x 







(57) 



Then for general £ G C the asymptotic expansion for x — > oo is different for 
£ G (1, oo) (case I) and for £ E C \ (1, oo) (case II). In the first case it its 
given by 



cry(x) = 2k\x + 2k\ — 



(58) 



4.r 



[(4fc 2 - 1) sin2si(a;) + 2fci cos2si(x) - 2fci (4fc 2 + l)] 



— 2 fci (4fcj 2 + l) 2 sin2si(.T) - Mi + C(x~ 3 ) 
.t + 01 + 2fci In (x) , 
where the parameters fci , tfij depend on £ as 

1 



si(x) 



2vr 



In(C-l) 



(59) 
(60) 



J r 2 (ifcx + i/2) 



and Mi is a constant to be determined. In case II the asymptotic expansion 
is 



a v [x) 



x (cot sn(x) + 2fcn) 
2) 



(61) 



4 ( 4fc n 



1 



(sinsn(x))^ 



2fcn cot sn(cc) 



(1 + 3fc 2 j) 2 cot s n (x) 



(sins n (a;)) 2 
in 



Mr 



(1 + S/cj^cot s n (x)+ 



+ 0(x~ 2 ) 



(sin s n (a;)) 2 
sn(x) = x + 4>u + 2ku In (x) , 

where in, Mu are unknown constants. The parameters are no 

i 



2tt 



In (l-£) , 

r 2 ^„) 



(62) 
(63) 



1 There is a typo in Eq. (1.29c) of Ref. [JTj. The square bracket has to read 
[e+^-n 2 -!]. 
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Only for A = 0, £ € (l,oo) and case I applies. As mentioned above for 
A = the extra particle's Green's function becomes identical with the single 
particle Green's function of a system of identical hardcore Bosons. For that 
case Vaidya and Tracy pmi52"] derived an asymptotic expansion, whose first 
terms read 



y/x\ Sx 2 



cos(2x) 

4 



0{x- A ) \ . (64) 



This result can not completely be derived from Eqs. (|58p and ([5^]). The 
constant Goo « 7ry / e2 _1 / 3 A" 6 , where A « 1.2842 is Glaisher's constant, 
cannot be obtained from the asymptotic expansion of ay. Neither can the 
constant term in the square bracket. Rather Eq. fixes the constant M\ = 
2" 4 . 

For | A | > case II applies. We write fcn = lot, where 

a = - - arctan(A) (65) 

2 7T 

is a number between —1/2 and 1/2. The cases a = and a/0 have to be 
treated separately. By expanding Eq. (|62|) for a = up to leading order we 
obtain ay ~ a; cot (a;) — 1. This together with K = yields 

Gi(x) = ^^i^M _ (66) 

x 

This is the well known Green's function of a system of free Fermions. Likewise, 
it could have been derived from Eq. (|43|) in the limit A — s- oo. It is easily 
checked that a: cot (a;) — 1 is an exact solution to the Painlcve V differential 
equation (|54"|) to the boundary condition lim^-j.o °v = —x 2 /3. 
For a^O the asymptotic expansion yields 

+C(a;- 4 l Q l- min ( 4 l Q l' 1 ))) . (67) 

The exponent of the leading order term can be written as a function of the 
extra particle's momentum K. It interpolates between —1/2 for the extra 
particles momentum in the core of the Fermi-sea and —1 for the extra par- 
ticle's momentum right at the edge of the Fermi-sea. As a approaches zero 
the approximation in Eq. (|67[) becomes poorer and poorer. Finally for a — 
infinitely many terms contribute to the same order and sum up to the simple 
result ([6"6")l . Within the present approach the constants and M\\ of the 
asymptotic expansion can not be determined. The first subleading term in 
the expansion is of order min(4|a|, 1). 

In Fig.[5]Gi(?/ - ) is plotted for different values of A. The curves are gener- 
ated from the representation of G\{y~) as a Toplitz determinant (see Eq. ([97)1 
in the appendix [A"|l . The number of sea particles is N = 29. The expansion 
for large distances in the thermodynamic limit as given in Eqs. (I55I) . (I64|1 
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O.Ott 0.5tt I.Ott 1.5* 2.0tt 2.5jt 3.0n 3.5n 4.0tt 

x = k P y~ 



Fig. 5 Real part of Gi(y~) for N = 29. The values of A are A = 0, A « 1.1 and 
A = oo. The thin dot-dashed lines correspond to the large x expansion (black). The 
thick dotted line (black) corresponds to the small x expansion for A = 0. The full 
lines (online red) are obtained from Eq. (|97l) . 

and (|57| are compared with the finite N result. The constants G^ 1 ^ « 0.6 
and M\\ ~ were numerically approximated. In particular for the non — 
oscillatory curve for A = the asymptotic expression works remarkably well 
almost everywhere. 

3.2 Finite Interaction Strength 

Up to now we were unable to evaluate the Green's function in the thermo- 
dynamic limit for finite interaction strength. Nevertheless the representation 
PTj) of Gi(y~) allows us to deduce the small distances behavior. We expand 
Gi(y~) for small y~ . As starting point we take the quantities g nm defined in 
Eq. (|42[) . Expanding them up to the quadratic term yields 

gnm = (b n S nm - 2c) + (i(-k n b n 6 nm + 2c{k n + k m - A)) - 2c 2 ) y~ (68) 

+ — [— 2k^ l b n S nm + 4c(fc 2 + k~ n + k n k m ) 

-Ac(k n + k m )(A - ic)} (y-) 2 + O {{y-f) , 
where we have introduced the notation 



b n = L[{k n -A) 2 + c 2 ]+2c 



(69) 
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For the linear combination of the quantities g nm as it appears in the Eq. (|41| 
we obtain 

Qnra Qn+lm Qnm+X Qn+lra+1 

(70) 

= \{b n + b n+ i)S nm — b n 5 nm +l — & n +ll5ri+lm] 

- i [(b n k n + b n +i ik n+ i5 n+ im] y 

- 2 ( + &n+l^n+l)^nm ~ i>n^nra+l ~ ^71+ 1 + 1 <5 rl + 1 m . ] 

-2c(fc„ - fc„ +1 )(fc m - fe m+ i)) (2T) 2 + 0((j/~) 3 ) ■ 

We plug the expansion ([70)1 into Eq. (|41[) and write Gj(y~) formally in a 
power series up to the quadratic term 

Gx(y-) = 1 + iG^y- G[ 2) (y-) 2 + O {(y-f) . (71) 

Our numerical results (see Figs. [3] and H]) together with the results obtained 
in the hardcore limit indicate that Gj(y~) is nodclcss for K < fcp. Therefore 
the short distances decay of the Greens function is dominated by the length 
[Gj ] -1 / 2 to which we refer as correlation length. In Eq. (f7T|) we used that 
Gi(0) = 1 by definition. Comparing with Eq. (|4"Tj) we have 

(NIL |C W | 2 ) -1 = dct [(6„ + b n+1 )S nm - b n S nm+1 - b n+1 8 n+lm ] n m=1 ^ j V . 

'(72) 

The determinant on the r. h. s. can be evaluated yielding for the normaliza- 
tion constant 

/N+l \ 7V+1 

\C N \- 2 = LN\ (lib, (73) 

\j=l J j=l 3 

The evaluation of the first and second order coefficient is somewhat lengthy 
but straightforward. We obtain 



/N+l 1 \ 1 /N+l , T . , so N+l N+l 



By definition the antisymmetric part of G(y~) and Gj(y~) is purely imagi- 
nary while the symmetric part is always real. This is reflected in the expansion 
([71]) and in the expansion coefficients (|74| and (|75|) . In the thermodynamic 
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limit the coefficients fl74l) and d75l) translate to 



G W = « + fdk J K - k) 2 , (76) 

— A.'F 

/ +,£F 2 

G i 2) = 2v(k?,c,K) J dk (k-A)* + <? (7?) 
Vfe F 

+ fc F + &F 



2tt ,/ 7 [(k ~ A) 2 + c 2 }[(k> - A) 2 + c 2 ^ 

— kp — hp 

The integrals are elementary. However, the resulting expressions give little 
insight. Thus, wc do not state them here. Equating Gj and Gj 2 ' in the 
hardcore limit yields 

lim = k F K , lira G{ 2) = %(1 + 3K 2 ) . (78) 

Note that Eq. ([75)) is in perfect agreement with the results of the Painleve 
transcendental evaluation of G\{y~) in the hardcore limit (compare Eqs. (|55[) 
and (|56p ). On the other hand we obtain in the limit of vanishing interaction 
strength 

lim G T (1) = lim g\ 2) = (79) 

where yl| c= o = K has been used. According to the equation above we have 
Gi{y~) = 1 f° r c = 0. This is reasonable since then the distinguishable 
particle does not feel the presence of the Fermi-sea and is described by a 
perfect coherent state. Great simplification arises for the ground state where 
K = J a = A = 0. As obvious from Eq. ([75)1 the linear term then vanishes 
identically and the correlation length acquires the simple form 



[G[ 2) ]- 1/2 = -i- [#- ( — T7 — 7T- - 2Vl - carctanfc" 1 )) 

1 V2fc F L2tt \ v arctan(c- 1 ) ^ v v " 



-1/2 



.(80) 



In general the correlation length decays monotonically between the values 
given in Eq. (|79| and (|78|) as c varies from zero to 00. Comparing [g| 2 ^] -1 / 2 
for the ground state and exited states of the extra particle, it is found that 
it decreases faster with increasing K. 

In Fig. [6] wc show the interaction part of the Green's function for J a — 
K = and for an excited state with K = —0.6 for different values of the in- 
teraction strength. The particle number is N = 29. The decay is considerably 
faster for the excited state. 
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1 2 3 4 5 

x = fc F y _ 

Fig. 6 Real part of Gi(x) as function of x = fez??/ - for the ground state (solid 
lines) and for a state with finite K w —0.6 (dotted lines). The interaction strength 
c = c/A?f has the values c = (black), 6=1 (online red), c = 3 (online blue) and 
c — oo (online green). The plots are generated for N = 29 sea particles. 



4 Density-Density Correlation 

The equal time density-density correlation function for a particle of the 
Fermi-sea and the distinguishable particle has been calculated in Ref. [14] . 
However, there the results are restricted to the case J a = A = 0. Here we 
give a brief derivation of the corresponding quantity where the distinguish- 
able particle may occupy an exited state inside the Fermi-sea. 

The density-density correlation function is defined by the (N — l)-fold 
integral 



R(x,y) — N I dxi - I dx N -i \W(x\, . . . , xjv-i, x, y)\ 2 . (81) 



Using the wave function ([2]) the above expression is after some algebra rewrit- 
ten as 



N+1 N N+1 p t(k n -k m )z 

*M - E 1 2 - E A (Z)A . (z) + ° (X- 1 ) ( 82 ) 

n=l \A n {l) n?m An[Z) m[Z) 
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with z = x — y. In order to obtain Eq. (|82[) we used 

L 



l(kj - hi) 



^(kj-Af + c^t+OiN- 1 ). (83) 



Taking into account the normalization condition 

L L 

dx I dyR(x, y) = N 



(84) 



o o 



we obtain 



R{z) = 



L 



1 - 



V 



JV+l 

N £ A. (1)1 

n=l 1 



N+1 „i(fc„-fc m )z 



+ 0(A^ 1 ). (85) 



In the thermodynamic limit the summations can be replaced by integrations. 
Then by introducing the dimensionlcss quantities 



c = c/k F A = A/ki 



z = kpz 



R(z) = -R(z) 
P 



(86) 



the density-density correlation function acquires the form 

2 



R{z) = 1 



2«(1, c, K) 



+ i 



dk- 



i(k -A) + c 



OiN- 1 ). 



(87) 



We discuss Eq. (|87|) in three limits. First it is easily seen that in the limit of 
infinite interaction strength the density-density correlation function becomes 
identical with that of free Fermions i.e. 



lim R(z) = 1 



sin z 



(88) 



On the other hand the second term in Eq. (|87|) vanishes like c as c -> and 
hence it follows that 



lim R(z) = 1. 



(89) 



Finally we consider the behavior of the density-density correlation function 
for z — > 0. This is the density of sea-Fermions at the position of the dis- 
tinguishable particle. For hardcore Bosons the local correlations have been 
studied in Ref. [43]. According to Eqs. (|88p and (|89|) we expect a transition 
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from R(0) = 1 to R(0) = as c increases from zero to +oo. The integral in 
Eq. (|87[) can be evaluated and the resulting expression reads 



R(0) = 1 



2v{l,c, K) 



In 



l + i(l-A)/c 
+ A)/c j 



(90) 



For K = i.e. A = this simplifies to Eq. (|5T|) . In Fig. 2] we show the plot 
of R(z) as function of the distance z and the interaction strength c. 




While for small distances density-density correlation function decays mono- 
tonically from one to zero as c increases this behavior changes at large dis- 
tances and R(z) exhibits a local minimum as function of c. The corresponding 
plot for the ground state is qualitatively similar to that of Fig. @] Quanti- 
tatively they differ in the form of the "edge" which marks the transition 
between monotonic and non-monotonic behavior in c. Compared with the 
ground state this edge becomes more and more smeared out as K increases. 



5 Summary &: Conclusion 

For a particle interacting via a (^-potential with a Fermi-sea, we calculated 
expectation values as well as the extra particle's equal time Green's functions 
and the density-density correlation function in dependence of its momentum. 
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The many-body wave function could be expressed as a determinant. This 
allowed us to obtain a compact expression for the Green's function even for 
finite interaction strength. 

For infinite interaction strength the thermodynamic limit can be per- 
formed and the Green's function is a solution of Painleve V, with a short 
distance boundary condition, which depends on the particle's momentum. 
Using the connection formulae for Painleve V we found an algebraic decay 
for large distances as x with f3 = (K 2 + l)/2, where K = K/k F is the 
extra particle's momentum on the scale of the Fermi-momentum. The parti- 
cle undergoes a transition from a free Fermion (K = 1) to a hardcore Boson 
(K = 0) as its momentum moves from the edge to the core of the Fermi-sea. 

Up to now we were unable to perform the thermodynamic limit of the 
Green's function for finite interaction strength. However with the determi- 
nantal structure of the Green's function presented here, an analytical way to 
perform this limit even for finite interaction strength is not out of reach. 

Yet our numerical simulations together with the small distance analysis 
seem to indicate that the interaction strength does not change the Green's 
function quantitatively but rather renormalizes the parameters. We conjec- 
ture that the relation of the Green's function to solutions of Painleve equa- 
tions holds beyond the hardcore limit. Such a relation indeed exists for hard- 
core Bosons and was employed [4"4"PX] to calculate the first order corrections 
in 1/c to the hardcore result (see Eq. (jM)) ). 
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cial support by Studienstifutung des deutschen Volkes. HK acknowleges financial 
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A Representation as Toplitz Determinant 

For finite interaction strength the determinant in Eq. (|40p representing the Green's 
function is not a Toplitz determinant due to the non constant diagonal entries. 
However, in the limit c — > oo the Green's function has a rep rese ntation as Toplitz 
determinant. To reveal it we write the quantities gji in Eq. (|42|) as 

Em gji = Lc V l( ^ +ni+1+2arctan(A)/ ^ ( (A 2 + l)5 jt - 2(1 + tA) sin M \ (91) 

c->oo V TTUjl I 

Here we use the notations 

t=^(y-y')>0 , n jl =j-l (92) 

and furt hermore assume that the quantum number rij are given by the set in 
Eq. (|15|1 Now the diagonal entries are constant, i. e. indepe nde nt of kj. For the 
linear combination of the quantities g nm as it appears in Eq. (|41[) this yields 

lim ( 9j , - gn+i - g j+ u + g j+ u + i) = L c 2 e-^+ n ' +1+2 ^ t ^' )t ^ l (t, A){93) 
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where we have defined 

4 sm(nji\t\) cos(i) 



7ji (t,A) = (A 2 + l)(2cos(t)(^ -<5 3 -i+i -Sji-i) + 



"ji 



i -b\ i nf 81 ^^'- 1 )!*!) | sin (K' + 1 )l t D ^ (94) 

7r V n j; - 1 Uji + 1 ) ' 

The full Green's function (see Eq. ([41])) then reads 

lim G(t) <x e 21 det [ 7j , (t , A)] 1=1 w . (95) 

c — ^-oo J ' 1***1 

Employing the normalization condition G(0) — \jL this evaluates to 

1% arctanf \)t / ty 

}^ G{t) = L(^ + l)(A^ + l)^ det[7ji(f ' A) W....,iv- (96) 
Correspondingly the interaction part G\{t) reads 

— 2iN arctan(A)£/ tt 

}^o Gl ^ = (N + 1)(A" + 1)-^ d6t [7ji(t ' A)] ^'= 1 w • (97) 

For A = Eq. (|96[) is equivalent to the r epresen tation of the Green's function 
for hardcore Bosons as Toplitz-determinant |30U27| . If on the other hand A — > oo 
we obtain from Eq. (|96l) 

it 

lim lim G(i) = det [(2 cos(t)^-, , - 5 jil+1 - !=1 jv . (98) 

£ N (t) 

Expanding the determinant in the equation above yields the following recursion 

C N (t) = 2cos{t)C N - 1 {t) - C N -2(t). (99) 

With help of the relation 2 cos(x) sin (a;) = sin(x — y) + sin(x + y) we see that its 
solution is given by 

= sin((iV + l)^ 
Sin(t) 

Consequently Eq. (|98|) acquires the form 

lim lim G{t) = f sin((jV + l) t ) ^ 
a-s-ooc^oo w L(JV + 1) sin(t) v ' 

Analogously one obtains for the interaction part 



lim lim Gi(i) = — u . . . ' ' . (102 



Eqs. (|101|l and (|102|l correspond to the free Fermion result. 
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